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Abstract. We study the classical correlation (CC) and quantum discord (QD) 
between two spin subgroups of the Lipkin-Meshkov-Glick (LMG) model in both binary 
and trinary decompositions of spins. In the case of bipartition, we find that the classical 
correlations and all the quantum correlations including the QD, the entanglement of 
formation (EoF) and the logarithmic negativity (LN) are divergent in the same singular 
behavior at the critical point of the LMG model. In the case of tripartition, however, 
the classical correlation is still divergent but all the quantum correlation measures 
remain finite at the critical point. The present result shows that the classical correlation 
is very robust but the quantum correlation is much frangible to the environment 
disturbance. The present result may also lead to the conjecture that the classical 
correlation is responsible for the singularity behavior of physics quantities at critical 
points of a many-body quantum system. 
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1. Introduction 

In a many-body quantum system, the interplay of various energies leads to different 
phases. When one of the energies becomes dominant over all the others by varying either 
adjustable interaction constants or applied external fields in the Hamiltonian, the system 
undergoes a phase transition and some observable display the singular behavior at a 
critical point. Since this phase transition occurs at zero temperature and is induced 
purely by quantum fluctuations in the system, in contrast to usual phase transitions 
induced by the thermal fluctuation, it is called quantum phase transition (QPT) [TJ. 
During the last decade, the QPT has attracted a lot of attention and become an 
important research domain pQ. On the other hand, quantum many-body systems have 
genuinely "quantum" correlations or entanglement in contrast to classical correlations 
[21 [3]. Therefore, it becomes natural to connect QPTs to quantum entanglement. At 
present, many measures for entanglement have been proposed such as the relative 
entropy [3], the concurrence (5], the entanglement of formation [6], the logarithmic 
negativity [7] and so on. As observables for identifying QPTs, those quantities indeed 
display the singular behaviour at zero temperatures [H EJ [10]. When calculating the 
measures of entanglement, one need to divide a system into several subsystems and 
then investigate quantum correlations between the subsystems. For a bipartition, two 
subsystems are of complementary parts one another and the entire system is always in 
a pure state. For a more multi partition, however, any two subsystems no longer forms 
a whole system and are in general in a mixed state. In this case, the other parts play 
a role of environments to the two subsystems under consideration and may strongly 
affect the critical behavior of quantum correlations between the chosen subsystems. In 
fact, Osborne et al. [UJ investigated the two-spin entanglement in the XY spin model 
and found that the entanglement remains finite and displays a peak at the critical 
point. Vidal et al. [12] studied the entanglement of a L spins block and found that the 
entanglement entropy displays a logarithmic divergence for large L. Recently, Werlang 
et al. [13] showed that the quantum discord (QD) [HI [15] and the entanglement of 
formation (EoF) between nearest-neighbors spins in an XXZ infinite spin chain at finite 
temperatures are no longer divergent but become finite at the critical point. 

The Lipkin-Meshkov-Glick model (LMG) is one of few solvable many-body systems. 
In recent years, a lot of efforts have been devoted to the study of quantum correlations 
such as the entanglement entropy, the concurrence and the logarithmic negativity in the 
LMG model [IHl HZl HH [191 [2QI EH [221 [23l [Ml I25l [26]. Latorre et al. [18] investigated 
the entanglement entropy in the LMG model and found that the entropy displays a 
singularity at the critical point. Morrison et al. [21] studied the dynamical QPTs with 
the spin-spin entanglement in a dissipative LMG. Orus et al. [23] investigated the 
many-body entanglement in the LMG model and showed that the critical scaling laws 
for the single-copy entanglement and the global geometric entanglement are equivalent. 
Wichterich et al. [27J found that the logarithmic negativity of two macroscopic sets 
of spins becomes finite at the critical point in any tripartition of mutually interacting 
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spins described by the LMG model whereas it displays a logarithmic divergence in a 
complementary bipartition. 

From previous studies, a question is raised which correlations are responsible for the 
divergent behavior of many-body systems at critical points. To clarify this question, one 
need to distinguish classical correlation from quantum one since the later one is much 
sensitive to the disturbance of environments. By recognizing the discrepancy between 
quantum extensions of two equivalent expressions for the classical mutual information 
[28] . Olliver and Zurek [TJ1 EE] introduced quantum discord (QD) that is sufficiently to 
qualify the total amount of quantum correlation including entanglement in a composite 
system and classical correlation (CC). In this paper, in order to analytically and clearly 
answer the raised question, taking the LMG model as an example, we compute the QD 
and the CC of two macroscopic sets of the mutually interacting spins in the cases of 
both bipartition and tripartition. For completeness and comparison, the EoF and the 
logarithmic negativity (LN) are also computed. We find that at the critical point both 
the CC and quantum correlations, including QD, EoF and LN, are always divergent in 
the bipartition. However, in the tripartition the quantum correlations remain finite and 
the CC is still divergent. The result may lead to a conjecture that the singular behaviour 
of observables of a many-body system with finite temperatures at critical points comes 
from the CC divergency. 

This paper is organized as follows. In section 2, the model is introduced. In section 
3, correlations in a bipartition setting are studied and detailed discussions are given. In 
section 4, correlations in a tripartition setting are investigated. Finally, a brief summary 
is given in section 5. 



2. The Model 

The LMG model describes a collection of mutually interacting N spins-1/2 on the 
x-y plane with an external field applied along the z direction. The Hamiltonian of LMG 
model reads 

1 N N 

h = ~m (*M+7*W) (i) 
i«j)=i i=i 

where of ((3 = x, y, z) are the Pauli matrices for a spin at position i, N is the total 
number of spins, < 7 < 1 is an anisotropy parameter and h is an external magnetic 
field applied along the z direction. 

In terms of the total spin operators Sp = J^i^i/^, the Hamiltonian ([[]) can be 
rewritten as 

H = -^(sl + 1 Sl)-hS z . (2) 

The ground-state properties of the LMG model have been found by use of a mean-field 
approach [191 [29]. The LMG model undergoes a second-order phase transition at h — 1. 
For h > 1, the ground state is a symmetrical and fully polarized state where all the 
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spins are along the external field direction. For h < 1, the corresponding ground state 
is two-fold degenerate [HI [29] . 



3. The quantum and classical correlations in a bipartition 

In this section, we divide the N spins into two groups and investigate the ground- 
state quantum and classical correlations between the two spin groups. To do so, we 
first need to determine the lowest energy state of the LMG model for a given external 
field. Thus, a rotation transformation to the total spin operators around the y axis is 
introduced as follows 



( s x \ 
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Sy = 




\s z ) 


\ 



COS^o 





sin#o \ 

1 



(S x \ 



(3) 
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\S Z J 

In ([3]), 9q stands for the value of the rotation angle which is chosen to make the 
expectation of the Hamiltonian (2) in the ground state (H) be minimum. 6 = 
for the symmetrical phase with h > 1, and 8q = arccos/i for the broken phase with 
< h< 1 pH[29]. Substituting Eq. (\3§ into P, one obtains 



H = - — cos 6»o sin 6 (S + S z + S„S Z + S Z S + + S Z SJ) 

~ Tm (( cos2 00 - 1) ( §2 + + §2 -) + ( cos2 6 ° + 7) (s + s- + S-S+)) 



2N 
1 

AN 

— sin 6 S; + — - — ( S + + S_ 



h cos 6 S Z . (4) 
we have used the relations S x = (^S + + S-j jl and 



When working out Eq. 

S v =(S+-S-)/(2i). 

We now split the N spins into two groups and consequently write the total spin 
operators as Sp = Sp + . In the Holstein-Primakoff representation [19], the spin 
operators Sp k \k = 1, 2) for each of the spin groups can be written as 
gW = Nk /2 - 4 



(fc) 



S 



S { _ k) 



N, 



<X k <X k , 

t A 1/2 
a' k a k J a k , 

sl/2 



(5) 
(6) 

a k [N k - a k a k J , (7) 

where a k and a k are bosonic annihilation and create operators and N k denotes the spin 
number in the kth group under the condition N = Ni + N 2 . 

Upon substituting Eqs. ©-([7]) into (j4]) and expanding if as a series of powers 
1/N k , and keeping the lowest order, one obtains 



H= - 



m 



7 



AN 

m 2 + 7 



4N 



+ 1 



N iC l\ + y/N^a^ + ^N^^ax + N 2 a 2 2 + h.c. 

(Ni^al + ^NiNzaial + \J ' N 1 N 2 a 2 a\ + N 2 a 2 a\ + h.c. 
+ /imj (a\ai + a 2 a 2 ^j + Const, 



(8) 
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with m = cos#o- The Hamiltonian (jSJ) is of a quadratic form of the bosonic annihilation 
and creation operators. It can be diagonalized by introducing the Bogoliubov 
transformation 



a 2 



cosh ® 61 + sinh ®b\ ) + b 2y ^, 



^cosh 



2 h + sinh — b[ ) ^fr 2 - fay/n, 



(9) 



(10) 



where bi(i = 1, 2) are new bosonic operators and Tk = Nk/N with J2k T k = 1- If choosing 

tanh0 = -s/r, (11) 

with 5 = 7 — m 2 , and r = 2/im — 3m 2 + 2 — 7, the Hamiltonian (jSJ) can be written in 
the diagonal form except an irrelevant constant to the present investigation 

H = A 1 b\b 1 + A 2 blb 2 , 



where 



Ai = — (2hm — 3m 2 + 2 — 7^ cosh B — (m 2 — 7) sinh G 
A 2 = - (2hm - 3m 2 + 2 - 7 



(12) 



(13) 



(14) 

The Hamiltonian (T121) represents two independent harmonic oscillators which 
ground state \ipo) is a Gaussian state, defined as 6, j^o) = 0. Therefore, the ground 
state of the LMG model can be fully characterized by the covariance matrix with the 
elements = (ipo\ \ipa), where R = (xi,pi,X2,p 2 ) with canonical coordinates 

Xk = (a\ + OfeJ / V% and momenta pk = i (a\ — a^j / y/2. By use of the Bogoliubov 
transformation ([9]) and ( fTUl) . one can obtain the explicit expression for the covariance 
matrix 

/ Am + 1 A X ^F 2 \ 
BiTi + 1 5 lv ^ 
^lv 7 ^ Axt 2 + 1 

fiiv^ B x t 2 + 1 / 



V 



C\ G 2 



(15) 



where 



Ax = ^J(r-s)/{r + s) 



B x = yj(r + s)/(r-s) 



(16) 
(17) 



and Gi, C{ are 2x2 matrices. For the simplicity of expressions in the following, we set 
\J (r + s) / (r — s) = a. In the symmetrical (h > 1) and broken (0 < h < 1) phases the 
parameter a reads 

^ - 1) / (h - 7), h>l 

;i-/ i 2 )/(i-7) ! < < 1 



a 
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By performing a like-Bogoliubov transformation [30] , the covariance matrix (j!5p can be 
written in the standard form 

(a ci \ 
a c 2 



ci b 

V o c 2 o 6 y 

where a, b, c\ and c 2 are determined by 
a 2 = det Gi = A, 

b 2 = det G 2 = B, 

c\c<i = det C\ = C, 



with 



(19) 



ab — cA (ab — c 2 J = det T = D 



A = B = or 1 [an + (l - n)] (n + a (l - n)) 



c 



2 — a — a 1 



[l-nln, 



(20) 
(21) 
(22) 
(23) 

(24) 
(25) 



and D = 1. 

Correspondingly, the simplectic eigenvalues of the covariance matrix ( fl~9l) are given 

by 



(M ± VM 2 - 4D) /2, 



(26) 



with M = A + 5 + 2C. 

According to Ref. |31j . the classical and quantum correlations of the Gaussian 
state characterized by the covariance matrix (fT9~j) can be respectively calculated by the 
formula 



cc = f (Va) - f (v 7 ^ 

QD = f (VB) - f - / + / (v^) , 



(27) 
(28) 



where E min = (2C 2 + {B - 1) (D - A) + 2 |C| ^/C 2 + (5 - 1) (D - A)) j [B - l) 2 for 
(D-AB) 2 < (1 + B) C 2 (A + £)), and when it comes to other cases E min 



is determined by \ AB - C 2 + D - ^/C 4 + (D — AB) 2 - 2C 2 (AB + D)j / (2B), and 

/(*) = (*?) In(i?)-(s=i) ln(==i). 

Upon substituting Eqs. ( p4^) - (!26|) into Eqs. (|27|) and (!28l) . the explicit expressions 
for the CC and the QD can be obtained as 



(29) 



2 2 2 2 

It is noted that the expression of the CC and the QD are the same as that of the 
entanglement entropy obtained in Ref. [TO] . 
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From fll9p . we can also obtain the entanglement of formation (EoF) between the 
two divided spin groups [6] 

EoF = f (A) , (30) 



where 



/ (A) = c + (A) log 2 (c + (A)) - c (A) log 2 (c_ (A)) , (31) 

c±(A) = (A- 1/2 ±A 1/2 ) 2 /4, (32) 

A = a - c, (33) 

c = a-^VC-l + a^l-rOn. (34) 



From Eq. f fl9|) . one can also calculate the logarithmic negativity (LN) [7J which 
has been obtained by Wichterich et al in Ref . [27] . 

In Fig. 1, the various correlations such as CC, QD, EoF and LN are plotted against 
the external field h. It is clearly shown that all the correlations between the two spin 
groups diverge at the critical point h = 1. By comparing the two figures, one may find 
that the anisotropic parameter 7 has little impact on the singular behaviour of these 
correlations and the CC, QD, EoF and LN display the same divergency at the critical 
point although they describe the different correlations. In Fig. 2, the correlations versus 
the external field are shown for different divisions of bipartition. It is observed that the 
divergent behavior of the correlations at the critical point is hardly affected by the 
bipartition way. 

In order to analytically investigate the critical behavior of the CC, the QD and the 
EoF, we expand Eqs. ( 129]) and (130]) at the critical point (h = 1) and obtain 

CC = QD = -- In (h - 1) + - In (1 - 7) + - In n (1 - n) - In 2. (35) 

and 

EoF = -1 log 2 (h - 1) + ~ log 2 (1 - 7) + i log 2 7i (1 - 7i) - 1. (36) 

Eq. (135]) shows that when reaching the critical point the CC and QD diverge as 
— |ln(/i— 1) which is consistent with that appears in Figs. 1 and 2. Interestingly, 
the singular behavior of the CC and QD is really the same as that of the logarithmic 
negativity [27], the entanglement entropy [10] and the single-copy entanglement [23J. 
From Eq. (I3"61) we know that the EoF diverges as — |log 2 (h — 1) at the critical point 
and behaves slightly different from the CC and the QD. 

Based on the scaling hypothesis proposed in Refs. [101 HH] , the finite-size scaling 
behavior of the CC and QD can be straightforwardly extracted from Eq. (|35|) 



CC = QD~ IlniV+ Jin (l- 7 ) + Jinn (1- 7i). (37) 
00 2 

This finite-size scaling behavior is identical to that of the logarithmic negativity 
[27], the entanglement entropy [TO], the geometric entanglement and the single-copy 
entanglement [23]. Therefore, all the correlations between artificial divided two parts of 
the mutually interacting spins in the LMG model obey the same critical scaling law. 
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4. The classical and quantum correlations in a tripartition 

In this section, we divide the mutually interacting N spins in the LMG model into 
three groups, each of which has iVj spins under the condition N = Ni + N 2 + N 3 , and 
investigate correlations between any two groups. In this case, if we consider the first 
and third groups, we need to trace the variable of the second group. Thus, the spins 
in the second group plays a role of the environment to the spins in the first and third 
groups, and the spins in the groups under consideration is generally in a mixed state. 

Following the same procedure as shown in the preceding section, we can diagonalize 
the Hamiltonian ([1]) and obtain the ground state of the LMG model, from which the 
density matrix of the ground state can be built. By tracing the density matrix over the 
variable of spins in the second group, one can obtain the reduced density matrix for 
spins in the first and third groups. It is obvious that the reduced density matrix is also 
of a Gaussian state. The covariance matrix of the reduced density matrix is found to 
be [27] 



( A\Ti + 1 A ly ^ \ 

Bm + 1 o B x ^n 

Av^3 o M + i o 

V B X ^T Z Bm + 1 / 



(3? 



If one sets T\ = t% = r < 1/2, the standard form of (1381) is the same as (fl9j) which 
elements are determined by 

A = B = a' 1 (or + (1 - r)) (r + a (1 - r)) , (39) 

C = -of 1 (a - 1)V, (40) 

D = a~ 1 (a + 2(a-l) 2 (t-2t 2 )) (41) 

according to Eqs. (20)-(23). In this case, the reduced density matrix is of a symmetrical 
Gaussian state [30J . The symplectic eigenvalues of Eq. (j3"81) are found to be 

v- = 1, 

i/+ = a- 1 ' 2 yja + 2 (a - l) 2 (r - 2r 2 ) 1 ^ ' 

Upon substituting Eqs. p)]) -(|4"2 ]) into Eqs. (I27jl - (j25]l . one can work out the the CC 
and the QD between spins in the first and third groups with E mm = 1 for h = y/j, and 

E min = (^_ 2a 2 _ Adar _ d ^ + ( a _ gj a ) r 2 + 2rf 2 r 3 + | rf 3/2| ^ + y T 2 ^ _ ^ j ^ for 

other circumstance with d = (a — l) 2 and // = 2a (a (r — 1) — r) (1 + (a — 1) r). Since 
the analytical expressions for the CC and the QD are much lengthy, we here have to 
give up to explicitly write them out. 

The entanglement of formation (EoF) can be obtained from Eq. f)3Up with 

A = ^(y/A-h) (VA~k 2 ), (43) 



9 



where 



fci 



( a _ I) 2 r2(l + (a-l)r2 (44) 



\ a (a + r — ar) 



fa-lVV 



k 2 = ± • (45) 

a * k\ 

In Fig. 3, the various correlations for an equal tripartition T\ = t% = 1/3 are 
plotted as a function of the magnetic field h. It is clearly observed that the CC diverges 
whereas all the quantum correlation measures such as QD, EoF and LN remain finite 
at the critical point. 

To clearly look into the behavior of the CC and QD at the critical point, we expand 
the analytical expression of the QD at h — 1 and obtain 



QD = In v - ' + £ In I v v ' ' - I . (46) 




Eq. (T46|) shows that the QD indeed remains finite at the critical point. Moreover, the 
value of the QD is irrelative to the anisotropy parameter 7. This universal character is 
much similar to that of the logarithmic negativity as found in Ref. [27J. Eq. f|46|) also 
shows that when r approaches to 1/2, QD diverges as that in the bipartition setting. 

In the similar way, we can find the analytical expression for the classical correlation 
of the LMG model around the critical point 



CC 




(47) 



In contrast to the QD, the classical correlation between the two spin groups diverges as 
— \ In (h — 1) at the critical point. This divergent behavior is the same as that obtained 
in the bipartition setting. 

In quantum information theory, the total correlation of a bipartite quantum system 
is measured by the mutual information [321 Qualitatively, the total correlation 

equals to the QD plus the CC. From the present result, it is very clear that in a tripartite 
setting the classical correlation is responsible for the divergency of the total correlations 
at the critical point [27]. 

The critical behaviour of EoF can also be investigated from Eqs. (l3~0j) . (1431 - ( 1451) . 
However, the analytical expression of it is too lengthy to be explicitly written here. We 
just give the numerical results. In Fig. 4, the EoF and QD are plotted as a function of 
the partition parameter r at the critical point h = 1. It clearly shows that when r is 
less than 1/2 the EoF and QD remain finite. When r reaches 1/2 and the tripartition 
reduces to the bipartition, the EoF and QD go from finite to infinity. 
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5. Summary 

The Lipkin-Meshkov-Glick (LMG) model describes a collection of mutually 
interacting spins-1/2 in an external magnetic field. By dividing spins of the LMG 
model into two or three parts, we study the classical correlation (CC) and quantum 
correlation measures such as the quantum discord (QD), the entanglement of formation 
(EoF) and the logarithmic negativity (LN) between the two spin groups. In the case of 
bipartition, where the two spin groups are complementary and their ground state must 
be of a pure state, we find that the classical correlations and all the quantum correlations 
are divergent in the same singular behaviour at the critical point of the LMG model. 
In the case of tripartition, however, the classical correlation is still divergent but all the 
quantum correlation measures remain finite at the critical point. In a tripartition, the 
spin group traced out plays a role of environments and the other two spin groups are 
general in a mixed state. The present result shows that the classical correlation is very 
robust but the quantum correlation is much frangible to the environment disturbance. 
In the real situation, a many-body quantum system is unavoidably to be coupled to 
its surroundings and is in a mixed state. Therefore, the present result may lead to the 
conjecture that the classical correlation is responsible for the singularity behaviour of 
physics quantities at critical points of a many-body quantum system. 
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Figure Captions 

Fig. 1 Various correlations as a function of the magnetic field h for the bipartition 
with Ti = 1/3. The symbols shown in the inset of Fig. 1(a) are applicable to the curves 
of Fig. 1(b). 

Fig. 2 Various correlations as function of the magnetic field h for the different 
divisions of bipartition with T\ = 1/2, 1/6, 1/100 and 7 = 0.5. The symbols shown in 
the inset of Fig. 2(a) are applicable to the curves of Fig. 2(b) and Fig. 2(c). 

Fig. 3 The various correlations as function of the magnetic field h for an equal 
tripartition T\ = t 2 = r 3 = 1/3. The symbols shown in the inset of Fig. 3(a) are 
applicable to the curves of Fig. 3(b). 



Fig. 4 The EoF and QD versus the partition parameter r at the critical point 
h = l. 



